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Abstract 

The lytic-lysogenic switch of the temperate lactococcal phage TP901-1 is fundamentally different from that of 
phage lambda. In phage TP901-1, the lytic promoter Pl is repressed by CI whereas repression of the lysogenic 
promoter Pr requires the presence of both of the antagonistic regulator proteins, MOR and CI. We model the central 
part of the switch and compare the two cases for Pr repression: the one where the two regulators interact only on 
the DNA, and the other where the two regulators form a heteromer complex in the cytoplasm prior to DNA binding. 
The models are analyzed for bistability, and the predicted promoter repression folds are compared to experimental 
data. We conclude that the experimental data are best reproduced the latter case, where a heteromer complex forms 
in solution. We further find that CI sequestration by the formation of MOR:CI complexes in cytoplasm makes the 
genetic switch robust. 
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1. Introduction 

Phenotypic variability under homogeneous condi- 
tions can readily be obtained by interlinking multiple 
gene regulatory pathways. Several well-characterized 
examples of phenotypic variations are known to be im- 
portant for different developmental process of bacteria, 
such as the presence of persister cells in Staphylococ- 
cus aureus and E. coli, development of natural compe- 
tence and sporulation in Bacillus subtilis, and the choice 
between lytic or lysogenic growth of temperate bacte- 
riophages [1, 2]. Two distinguishable phenotypes may 
originate from a bistable system, i.e. a system that can 
toggle between two alternative stable steady-states ]3|]. 
Infection of bacteria by temperate bacteriophages pro- 
vides a classical example of the possibility to choose 
between two alternative modes of development. 

The bacteriophage lambda infecting Escherichia coli 
has been subjected to decades of intensive study, mak- 
ing the lytic-lysogenic switch one of the best understood 
gene regulatory systemsJUl^,!^]. The bistability of the 
lambda switch is obtained from a double negative feed- 
back mechanism, where two repressor proteins directly 
repress transcription of the other repressor gene. This 
system has a stable state with one promoter on and the 
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other off, and vice versa for the other stable state. Once 
either state has been established, it would persist indef- 
initely or until some trigger stimulus forces the system 
to switch to the other state. 

The genetic switch of the temperate lactococcal bac- 
teriophage TP901-1 infecting Lactococcus lactis subsp. 
cremoris provides a regulatory system diverse from the 
lambda genetic switch. A previous study has demon- 
strated that a DNA fragment obtained from the temper- 
ate lactococcal phage TP901-1 shows bistability when 
introduced into Lactococcus lactis. The cloned DNA 
fragment contains the two divergently oriented promot- 
ers, Pr and Pl, and the two promoter proximal genes 
cl and mor |7](Fig[Tk). A knockout mutation in the 
mor gene showed that CI ensures tight repression of the 
Pl promoter and partially repression of the Pr promoter 
whereas a knockout mutation in the cl gene results in 
open states of both Pr and Pl, showing that MOR by it- 
self does not exhibit repression of either promoter||8l|9|]. 
Two types of repression has been shown: i) MOR- 
independent repression, which is responsible for repres- 
sion of Pl. The Pl promoter is repressed by cooperative 
binding of CI to the three operator sites Or, Ol and Od, 
by the formation of a CTDNA loop structure, ii) MOR- 
dependent repression, which is responsible for repres- 
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Figure 1: The regulatory circuit of the genetic switch from phage 
TP901-1 may be described as a double negative regulatory network, 
(a) Wild type genetic switch with Ol, Or, and Od sites. The relative 
positions of the cl and mor genes and the divergently oriented pro- 
moters Pr and Pl are also shown. The cyan box represents a putative 
Om site, (b) Modified switch investigated in this paper. Only one 
CI operator site, Ol, is present (dark green box), (c) Immune state 
of the modified switch. CI represses transcription from Pl approxi- 
mately 1000-fold by binding to the Ol operator site with transcription 
from Pr being allowed, (d) Anti-immune state of the modified switch. 
Repression of Pr depends on both CI and MOR, which repress tran- 
scription from Pr approximately 100-fold with transcription from Pl 
being allowed. 



sion of Pr and only occurs in the presence of both MOR 
and CI. This repression is suggested to occur through 
MOR and CI binding at a putative Om operator site H . 
Hence, the bistability of the genetic switch from phage 
TP901-1 may be described as a mixed feedback loop, 
where both of the antagonistic repressor proteins are in- 
volved in Pr repression. It is still not clear how MOR 
and CI collectively repress transcription from Pr, and so 
far there is no direct experimental evidence for interac- 
tion between CI and MOR. 

In order to understand the mechanism of switching 
in phage TP901-1, we here study a modified version of 
the cloned wild-type TP901-1 switch(FigOJ)). This con- 
struct contains only one of the three CI operator sites, 
Ol, which gives tight repression of Pl thus still sustains 
the bistable behavior of the construct. In the immune 
state, the Pl promoter is repressed approximately 1,000- 
fold, but high expression from Pr is allowed due to the 
absence of Or, that autoregulates transcription from Pr 
in the case of the wild-type switch [9] (FigJTJ;). In the 
anti-immune state, Pr is repressed approximately 100- 
fold but high expression from Pl is allowedH (Fig fTJl) . 

Generally speaking, it is not easy to predict the be- 
havior of a bistable switch without quantitative analy- 
sis because a bistable switch is a dynamical and highly 
nonlinear system. In the present case of TP901-1, even 
the modified version of the switch (Fig[TJ)) could in- 
volve a number of mechanisms, such as cooperativ- 
ity binding via homo/hetero-dimerization, sequestration 
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Figure 2: Repression of Pl though binding of CI2 at Ol is employed 
in both models (blue box). Repression of Pr is supposed to occur by 
direct binding of CI and MOR to Om (Model A) or by formation of 
a MOR:CI complex prior to binding to Om (Model B). The statistical 
weight is given for the different states of the operator sites. In Model 
A, it is assumed that binding of CI alone to Om is so week that the 
operator is only occupied when binding is stabilized by MOR, or vice 
versus if MOR binds first. 



via heterodimerization, intertwined loops of negative 
and positive feed back via protein interactions, etc. In 
such a situation, the only way to obtain any reliable re- 
sults is to perform quantitative analysis on specific mod- 
els. By confronting numerical results with experimental 
data, we can restrict possible mechanisms with plausi- 
ble parameters 

In this paper, we construct mathematical models for 
this modified bistable system based upon statistical me- 
chanics, examine their behavior in the steady states nu- 
merically, and compare the obtained repression folds 
with experimental observations[9]. We assume that 
Pl is repressed by the CI dimer binding to the opera- 
tor Ol since the operator has two inverted repeated se- 
quences II 1 OTl . On the other hand, the MOR-dependent 



repression of Pr is assumed to be brought about by the 
MOR:CI:DNA complex formation at the putative Om 
operator. Since the amino sequence of MOR shows high 
similarity to the DNA-binding Helix-turn-Helix domain 
of the repressor protein encoded by Escherichia coli 
phage 434 HH MOR is likely to be a DNA binding pro- 
tein, but the fact that MOR alone does not repress tran- 
scription from either Pl or Pr suggests that the DNA- 
binding affinity of MOR is negligible. CI alone does 
not repress Pr in this system due to absence of Or opera- 
tor. Based upon these observations, we test two different 
scenarios for Pr repression(Fig|2]i: i) repression through 
direct binding of CI and MOR to Om (Model A), and ii) 
repression through binding of MOR:CI complex to Om 
with the complex being formed in the cytoplasm prior 
to the binding to DNA (Model B). 

The MOR:CI:DNA complex formation in Model A 
may be regarded as an extreme case of Model B where 
the MOR:CI complex formation in cytoplasm is so 
weak that the complex is stabilized only when it binds to 
DNA, thus there is no substantial presence of MOR:CI 
complex in the cytoplasm. However, distinguishing 
Model A from Model B helps us to recognize two dis- 
tinct aspects in the Pr repression by the MOR:CI:DNA 
complex, i.e., co-operativity and sequestration; The lat- 
ter has been studied in silico as a possible mechanism 
for a genetic switch jl2l 1 1311 and demonstrated to pro- 
vide strong nonlinearity[14j |l5|]. In fact, our Model 
B is a reminiscent of one of the bistable switches ob- 
tained by the simulated evolution (Fig.3A in iS). We 
will demonstrate that CI sequestration by the MOR:CI 
formation in cytoplasm can make a robust bistable sys- 
tem and is actually a plausible switching mechanism for 
TP901-1. 



Theory 

The regulatory circuit in the present system consists 
of the two promoters Pl and Pr, which produce MOR 
and CI, respectively(Fig[lJ>). The promoter Pl is re- 
pressed by CI binding at Ol, thus the Pl activity is given 
by a function of the CI concentration as 



P L([C\]) = pL ■ / 0L ([CI]), 



(1) 



where pLo is the bare activity of the promoter Pl. The 
function /ol([CI]) represents the repression factor. In 
the absence of CI, there is no repression: /ol(0) = 1. 
The Pr activity, on the other hand, depends on the con- 
centrations of both CI and MOR. Accordingly, the Pr 
activity can be written as 

pfl([MOR], [CI]) = pR • / OM ([MOR], [CI]) (2) 



with pRo being the bare activity. The function 
/om([MOR], [CI]) is the repression factor due to 
the binding of MOR and CI at Om, and satisfies 
/ OM (0, [CI]) = / OM ([MOR],0) = 1. 

The promoter Pl produces MOR, and Pr produces CI, 
thus in the modeled feedback system, the total concen- 
tration for each protein, [MOR] to tai and [CI] to tai, is gov- 
erned by the dynamics equations, 

-r [MOR] total = — ( P L([CI]) - [MOR] tot A (3) 
at tm\ i 

^[CI] tot ai = ^(^([MOR], [CI]) - [CI] tota i), (4) 



where tm and re are the degradation times for MOR 
and CI, respectively. To simplify the notation, we have 
rescaled the promoter activities, Eqs.O and (0, by the 
degradation times, i.e. the promoter activities are now 
measured in terms of the steady state protein concentra- 
tions. 

In steady states, the production and the degradation of 
each protein should balance, therefore, the promoter ac- 
tivities and the concentrations of the expressed proteins 
in the cytoplasm should satisfy the steady state condi- 
tion, 



P L([Cl]) = [MORW, 
/?fl([MOR],[CI]) = [CIW 



(5) 
(6) 



Not all steady states are stable against small pertur- 
bations. A steady state is unstable if a perturbation 
drives the system out of the state; The stability should 
be determined by the dynamics equations, Eqs.® and 
(See supplementary material). If there are two 
stable steady states, the system shows bistability. 

We assume that the repression factors, /ol and /om 
in Eqs.lQ]) and (0, are given by the statistical weights 
at equilibrium that the corresponding operators are not 
bound by the regulators. This approximation holds 
when the time that RNA polymerase (RNAP) needs to 
start elongation after binding to DNA is much shorter 
than the time scales of binding/unbinding of RNAP 
and repression factors to the promoter/operator sites 
The equilibrium statistical weights depend upon 
the repressor concentrations, and their dependence is 
characterized by the Hill coefficient and the affinities of 
the repressors to the operator sites Il7l[l8i[l9ll . 

For the MOR-independent repression of Pl, we sup- 
pose that Pl is repressed by CI dimer binding at Ol, and 



that the dimers are formed in the cytoplasm before bind- 
ing. Thus, within the above approximation for the re- 
pression factor, the Pl activity is given by 

" i<[CI1) = ^TT[C^- l7) 

where the affinity Kol represents the CI2 concentration 
at which Ol is occupied for 50% of the time. Since 
the dimer concentration [CI2] is related to the monomer 
concentration [CI] as 

[CI] 2 



[CI 2 ] 



Kcu 



(8) 



with the dissociation constant Kch, Eq.(|7]i may be writ- 
ten as 

1 



pL([Cl]) = pLo 



1 + [Cl] 2 /K? 



(9) 



OL 



with the effective affinity 

£ol = ^K cl2 K 0L (10) 

for CI concentration. In Figf3] the activity of Pl as a 
function of [CI] is plotted by a green line. 

As for the MOR-dependent repression ofPn, we will 
examine two models. In Model A, monomers of MOR 
and CI may bind cooperatively at Om, but we do not as- 
sume any MOR:CI complexes formed in cytoplasm be- 
fore binding to DNA. In Model B, on the other hand, CI 
and MOR may associate in cytoplasm before they bind 
at Om. For both models, Pr is repressed by the forma- 
tion of the MOR:CI:DNA complex at Om. The important 
point in Model B is that the formation of MOR: CI het- 
eromers competes with CI dimer formation by seques- 
tering CI monomers . 

Model A 

We first consider a MOR:CI:DNA complex contain- 
ing one MOR and one CI protein as illustrated in Fig|2] 
Then, we can approximate the total concentration of 
MOR unit by the MOR monomer concentration, 



[MOR] toml = [MOR]. 



(11) 



The activity of the Pr promoter is repressed from the 
bare activity pRo by the statistical weight that the oper- 
ator Om is not occupied by MOR and CI, 

1 



pfl([MOR],[CI]) = pfl 



1 + [MOR] [CI] /(K 0M ) 2 

The affinity Kom is the concentration V[MOR] • [CI] 
where Om is occupied by MOR and CI for 50% of the 
time. 



(12) 



The steady state is determined from the steady state 
condition Eqs.© and © by eliminating the MOR con- 
centrations. With the help of Eq. lfTTI) . we obtain 



pR(pL([Cl]), [CI]) = [CI] 



total; 



(13) 



which represents the balance between the production 
and the degradation of CI. This can be solved graphi- 
cally by plotting the both sides as a function of [CI], 



pR(pL([Cl]), [CI]) 



= pRo 



1 + 



pL()Ko L /(Ko M ) 2 



([Cl]/K 0L ) + (*ol/[CI]) 



[CI]total = [CI] + 2 



[CI] 



, (14) 



(15) 



Eq. lfPfl) represents the Pr activity in the system where 
MOR is provided by Pl but [CI] is controlled externally. 
Note that the relative strength of the bare promoters, 
pLo and pRo, does not affect the system behaviors, such 
as bistability or repression folds, because there is no di- 
rect interaction between MOR and CI in this model. 

Model B 

In this model, a MOR CI heterodimer is formed in 
solution before it binds to the putative Om site to repress 
Pr (Fig|2|i. The activity of the Pl promoter is again given 
by Eq.© but the Pr activity is 



p/J([MOR],[CI])=/7*o 



1 



1 + [MOR CI] I Kom 



,(16) 



where Kqm now represents the concentration of the 
MOR CI heterodimer at which Om is occupied for 50% 
of the time. The concentration of the MOR CI het- 
erodimer is given as 



[MOR CI] = 



[MOR] ■ [CI] 

•^MOR CI 



(17) 



with the dissociation constant A'morci for the het- 
erodimer. 

The formation of the heterodimers couples the 
monomer concentrations of CI and MOR through 

[CI] to tai = [CI]+2[CI 2 ] + [MORCI], (18) 
[MORJtotai = [MOR] + [MOR CI], (19) 

which leads to the competition between the CI2 forma- 
tion and the MOR CI formation. Note that [CI 2 ] is still 
given by Eq.©. 

The steady state is determined as in the case of Model 
A; We consider the Pr activity as a function of [CI] 



when MOR is provided by Pl. From Eqs.lfTTb and ( fl9l ), 
[MOR] is expressed in terms of [CI] and [MOR] to tai, 



[MOR] 



[MOR] tota i 



(20) 



1 + [CI]/^MOR:CI 

and then, in the steady state where Eq.© holds, 
[MOR] tota i is given by the Pl activity with Eq.0. Then 
the steady state condition Eq.© for Pr becomes 



pR 



1 + [CI] /A'mORiCI 



, [CI] =[CI] t , 



(21) 



which can be solved graphically with the explicit forms 
for the both sides: 
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pLdCH) 
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1 + 
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pLp Kql/{Kqm) 2 

([CI]/£ 0L ) + (K J[Cl]) ' 1 + [CI]/XmOR:CI 



1 



(22) 



rcn 2 

[Cl] to tal = [Cl]+2±— L +pIX[CI]) 

where the effective affinities are 



[CI]/^MOR:CI 
1 + [CIJ/A'MORiCI 



^(23) 



^om = V^mor:Ci ■ ^om and ^ol = -\jKci 2 Ko L . (24) 

Eq. (l2TI > represents the balance between the production 
and degradation of CI in the system where MOR is pro- 
vided by Pl. Note that Model B reduces to Model A in 
the limit of large Zl"mor:Ci with Kom being kept constant 
as has been discussed at the end of Introduction. 



2. Results 

We numerically examine the steady states for the two 
versions of the models we have constructed (Figj2j. 

Model A 

In our first model, we study the possibility for bista- 
bility in the system where Pr is repressed by binding a 
MOR monomer and a CI monomer to Om without direct 
interaction between MOR and CI in the cytoplasm. The 
binding affinities for each of the proteins alone at Om 
should be negligible because the Pr repression requires 
both of the proteins. Hence, the affinity Kom in Eq. (fT2l 
may be considered as the effective binding affinity Kom 
for MOR and CI with very weak MOR CI formation, or 
as the resulting binding affinities from CI:Om, MOR:Om, 
and the interaction between the bound proteins. 

Figure [3] shows a typical example of the Pl activity of 
Eq.® (green line), the Pr activity of Eq. (TT4b (solid red 
line), and the total concentration of CI, or degradation 
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Figure 3: Model A. Promoter activities, pL and pR, and [CI] tMa i 
as a function of free [CI]. Curves represent Pl activity of Eq.j5J 
(green line), Pr activity of Eq. )14t (solid red line), and [CI] tot ai of 
Eq.{T3)(dashed red line). The filled red circle at the intersection be- 
tween the [CI]t tai and Pr curve represents the steady state. 



rate, of Eq.dTsb (dashed red line) as a function of free 
[CI] in the logarithmic scale. One can see that Pl (green 
line) is fully active and produces a lot of MOR at low 
[CI], whereas its activity is monotonically decreasing 
with increasing [CI], due to Pl repression by CI2 binding 
at Ol; Pl is virtually shut down beyond [CI] ~ Kol = 
^K Ch ■ K OL = 6 • IO- 4 . 

On the other hand, the Pr activity (solid red line) 
shows a more complicated behavior, i.e., Pr is open 
both at very low and high [CI] concentrations but re- 
pressed at the intermediate concentration. One can un- 
derstand this behavior by noting that the solid red line, 
Eq. lfPfl) . represents the Pr activity in the system where 
MOR is expressed from Pl under CI control; At low 
[CI], there is plenty of MOR due to high Pl activity, 
while at high [CI], no MOR is present. The Pr activ- 
ity is repressed only at intermediate [CI] because both 
MOR and CI are necessary for its repression by means 
of the MOR:CI:DNA complex formation. 

In the steady state, the production and the degradation 
of each protein should balance. Since CI is expressed 
from Pr and the degradation of CI is assumed to be pro- 
portional to the total concentration of CI, the steady 
states are identified as pR(pL([Cl]), [CI]) = [CI] tota i, 
namely, Eq.([T3l. We thus find the steady states at the 
intersection points between the pR and [CI] tota ] curves 
represented by Eqs.lfPfli (solid red line) and ( [13] ) (red 
dashed line), respectively. 

In Fig|3] only one intersection point is observed be- 



tween the two curves, showing that there is only one 
steady state solution. This uniqueness of steady state 
holds true for any given value of the parameters, be- 
cause the activity of Pr given by Eq. (fT4b never in- 
creases faster than proportional to CI, whereas [CI] to tai 
by Eq.(fT5b always increases faster or proportional to 
[CI]. Therefore, bistability is never realized in Model 
A with the assumption of binding of one MOR and one 
CI for repression of Pr. 

Variant of Model A 

Bistability may be obtained in Model A if a larger 
number of proteins are allowed to form a complex 
structure at Om. Suppose m MOR monomers and c CI 
monomers bind at Om to form the MOR,„:CF:DNA 
complex that represses transcription from Pr, then for 
the expression for the Pr activity, Eq.([T2l> should be re- 
placed by 

/?fl([MOR],[CI]) 



F 1 + [MOR] m [CI] c /(^OM) m+e ' 
while all the other equations remain the same. Using 
this for the left hand side of the steady state equation 
(O with [MOR] oc 1/[CI] 2 from Eq.© for large [CI], 
one can see that the largest slope of the Pr activity as a 
function of [CI] is 2m — cm the logarithmic scale. Since 
the slope of the plot of [CI] to tai given by Eq.([T5ll is be- 
tween 1 and 2, the steady state equation ( fT3l l can have 
more than two solutions with Eq.(l25b when 2m - c > 2. 
In the case of m = 2, we could obtain multiple solutions 
with c — 1 or 2, i.e., two MOR monomers binding to- 
gether with one or two CI monomers at Om. Examples 
for Model A with (m, c) = (2, 1) and (2, 2) are shown in 
Fig. |4] In each example, the intersections between the 
solid red line (the Pr activity) and the dashed red line 
(the total CI concentration) represent steady state solu- 
tions. One can see there are three solutions for each case 
in Figg] 

Dynamical analysis shows that the steady state in the 
middle marked by an open red circle is unstable against 
small fluctuations, and the states at the ends marked by 
filled red circles are stable (See supplementary material 
for detail). Full analysis requires Eqs.Q and 10} , but 
the stability may be understood in the following way; 
For the steady state in the middle, if the CI monomer 
concentration increases by fluctuation from the steady 
value of [CI], the CI production from Pr will increase 
more than the increase in degradation given by [CI] to t a i, 
as is seen in Fig0] where the solid red line of the Pr ac- 
tivity goes above the dashed red line of [CI] to tai, upon in- 
creasing [CI] from the middle steady state. This means 
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Figure 4: Model A with (a) two MOR monomers and one CI 
monomer binding to Om, (m, c) = (2, 1), and (b) two MOR monomers 
and two CI monomers, (in, c) = (2, 2). Promoter activities, pL and 
pR, and [CI] to tal as function of free [CI]. Curves are Pl activity of 
Eq.{D (green line), Pr activity of Eq.^25) (solid red line), and [CI] tota i 
of Eq. )15K dashed red line). The filled red circles at the intersections 
represent the stable steady state, and the green circles indicate the 
Pl activity in the steady states. The open red circles mark the un- 
stable steady states. The repression folds />K(open)/pi?(closed) and 
pL(open)/pL(closed) are approximately 200 and 1000 respectively. 

that such a fluctuation causes further increases of [CI], 
consequently, the state is driven out of the steady state. 
On the other hand, the steady states at both ends rep- 
resent stable states. A fluctuation towards larger [CI] 
leads to insufficient CI production in comparison with 
the CI degradation, i.e. the solid red line goes under 
the dashed red line as [CI] increases. This brings the 
system back to the original state, therefore, they are sta- 
ble. Thus, the system has two stable steady states, which 
leads to bistabilityQ 

The promoter activities in the two stable states can 
be determined from the graphic representation in FiglU 
The Pr activity is read from the ordinates of the inter- 
section (filled red circles) and the Pl activity is read off 
from the corresponding [CI] values of the intersection 
points (green circles), which allows us to estimate the 
repression folds for Pl and Pr between the two stable 
states. The state at the right represents the immune state 
with open Pr and repressed Pl, while the one at the left 
represents the anti-immune state with open Pl and re- 
pressed Pr . 

The relative activities between the two states should 
be compared with the promoter activities obtained from 
the in vivo measurements H; Pl is repressed approxi- 
mately 1,000-fold in the immune state and Pr approx- 
imately 100-fold in the anti-immune state. To repro- 
duce these repression folds in Model A with (m, c) = 
(2,1) and (2,2), we test the three parameters, Kch, Kol, 
and Kqm, representing the dimerization constant of CI, 



' Note that this simplified analysis is based on the assumption that 
Eq.{5} always holds, i.e. the response of Pl is much faster than that of 
Pr. For general case, see supplementary material. 
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Figure 5: Distributions of possible parameters for Model A. The 
plotted parameter sets are those that give the repression folds 50 < 
pR(open)/pi?(closed) < 200 and pL(open)/pL(closed) > 500 out of 
10 s tested parameter sets, which are randomly picked in the logarith- 
mic scale over the range of A"ol e (10~ 5 , 10~ 2 ), Kqu 6 (10~ 4 , 10), 
andJsr C i 2 e (lO^ 3 , 10 2 ). 

the effective binding constant of Cl-monomer at Ol, and 
the binding constant of the MOR:CI complex at Om, re- 
spectively, by setting the criterion that the Pr repres- 
sion fold should be in the range from 50 to 200 ( 50 < 
pR(open) / pR(closed) < 200) and the Pl repression fold 
should at least be 500 (/?L(open)//?L(closed) > 500). 
Fig|5] shows the distributions of accepted values of pa- 
rameters out of randomly chosen values in the logarith- 
mic scale, ^ol and Aom are narrowly distributed while 
Kq\ 2 are larger than 10" 1 for (m,c) = (2, 1) and larger 
than 10~ 2 for (2,2) in the unit of CI concentration at full 
activity of Pr. One can also see that the accepted values 
for Kqi 2 are much larger than A^ol- This suggests that, in 
order for Model A to work, CI must exist as a monomer 
and act by cooperative binding to form CI2 at Ol when 
repressing Pl. 

Figure [6] shows the parameters that satisfy the crite- 
rion only for the Pr repression fold versus resulting Pl 
repression fold (left two columns for Model A and right 
two for Model B). The vertical green lines are drawn 
at the Pl repression fold 500, thus only the plots on 
the right side of the lines should be accepted by the 
repression fold criterion. From the plots for Kqi 2 for 
Model A, one can see that the relatively high values for 
Kqi 2 in this model comes from the requirement for the 
large Pl repression fold. This can be understood as fol- 
lows; In order to achieve large repression fold for Pl, 



the difference in [CI] for the two steady states should be 
large, which in turn requires smaller slope in the [CI] to tai 
curve, namely larger Kch because the slope in [CI] to tai 
changes from 1 to 2 around [CI] ~ Kch (PigHJ. 

Model B 

Now, we consider the possibility that Pr is repressed 
by a MOR:CI complex formed in cytoplasm before 
binding to DNA. Examples for Model B are shown in 
Fig|7] where the Pl activity of Eq.© (green line), the Pr 
activity of Eq. d22l) (solid red line), and the total density 
of CI of Eq.(l23]l (dashed red line) as a function of [CI] 
in the logarithmic scale; All of them have three steady 
states. 

Striking difference from the case of Model A is that 
the [CI] t otai (dashed red line) can be non-monotonic. 
Therefore, for a given [CI] to tai within a certain range, 
there exist three possible states with different [CI]. This 
suggests that the bistability could be obtained for the 
system with Pl even if Pr were not regulated, namely, 
even if Pr would produce CI at a fixed rate within the 
range. This bistability is due to the MOR:CI heteromer 
formation in cytoplasm; Pl produces MOR, which se- 
questers its own repressor, i.e. CI, by forming MOR:CI. 
For a given [CI] to tah the state at low [CI] is the state 
where most of CPs are incorporated in MOR:CI het- 
erodimers due to the MOR produced by Pl, while the 
state at high [CI] is the state where most of CI is in the 
dimer with Pl being repressed. The state in the mid- 
dle is unstable. Such bistability is, of course, not the 
bistability observed in the experiments, but one can see 
that this feature of behavior in [CI] to tai (dashed red line) 
makes it easier to have three intersections with the Pr 
activity curve (solid red line) than in the case of Model 
A. 

According to the stability criterion we discussed, the 
steady state at both ends are stable while the state in 
the middle is unstable even for the system where both 
Pr and Pl are regulated. Full analysis, however, shows 
there are some cases where the states at both ends can 
be unstable although the stability criterion is correct for 
most cases (See supplementary material). We analyse 
the bistability based upon the stability criterion, ignor- 
ing the small possibility that the states at both ends 
could be unstable. 




Figure 6: Distribution of possible parameters versus Pl repression fold, pL(open)/pL(closed) for Model A with (hi, c)=(2,l) and (2,2), and for 
Model B with (m, c)=(l,l) and (2,1). The plotted parameters are those that satisfy the Pr repression fold criterion 50 < />K(open)/pi?(closed) < 200 
out of randomly chosen parameters from the region Kol, Kom, Kci 2 £ [10~ 7 , 1] for Model A and ATol. Kom, Kci 2 , ^MOR:Cl £ [10~ 7 , 1] for Model 
B. Vertical green lines are drawn at Pl repression fold 500, thus only the parameters that are in the right side of the lines are consistent with 
the experimentally obtained repression folds. The number of tested data are 10 6 for Model A with (m,c) = (2, 1), 5 X 10 6 for Model A with 
(m, c) = (2, 2), and 2 X 10 s for Model B. Note that the effective affinity £ M for Model B is defined by K'^ = K 0M K'^ R l a . 
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(m,c)=(l, 1), pL =l, P R =l 
K OL =3.92xl0 4 . K a =5.00X10 
K OM =5.74x!0' 4 , ^MORCf- * "' 



(m, ch(l, 1), pL =l. pR„=0.1 
K OL =3.21xW s . K a =6.45X10' 
K OM =2.18 xia 4 , K mma =3.20 x/0' 2 




(m, c)=(2. 1). pL„=l. pR =l 
K OL =7.69xia 4 . K cl =1.55x10-' 
K OM =4 .18xlff 3 . K MOR . c ,=1.90xlO- 2 



(m,c)=p,l), pL a =lpR =0.1 
K OL =3M>xl0 4 . K a =5.39X10-' 
K tm =4.50 xlO 1 . K mm a =6.57 xW 2 




Figure 7: Model B with MOR CI binding to Om, (m,c) = (1, 1). 
The bare promoter activities, pLg and pRg, are equal in left panel 
while the bare activity pLo is 10 times stronger than pRo in the right 
panel. The upper graphs show examples of the promoter activities, 
pL and pR, and [CI] tota i as a function of [CI]. The repression folds 
of Pl and Pr are approximately 1000 and 50, respectively in the left 
figure, and approximately 60,000 and 100 in the right figure. The 
middle graphs show the ratios of MOR units in the forms of monomer 
and heteromer, and the lower graphs show the ratios of CI units in the 
form of monomer, dimer, and heteromer. One can see that MOR:CI 
and CI2 compete for the free CI units in the intermediate concentration 
range of [CI]. 



Figure 8: Model B with MOR 2 CI binding to Om, (m,c) = (2, 1). 
The bare promoter activities, pLo and pRg, are equal in left panel 
while the bare activity pLo is 10 times stronger than pRo in the right 
panel. The upper graphs show examples of the promoter activities, 
pL and pR, and [CI] tota i as a function of [CI]. The repression folds 
of Pl and Pr are approximately 1,000 and 50, respectively, in left 
figure, and approximately 1,000 and 60 in right figure. In contrast 
to the case with (m,c) = (1, 1) in Fig0 the heteromer MOR2:CI 
is not formed at high [CI] because the Pl promoter is closed faster 
than the MOR2:CI is formed. Note that the effective affinity of Om is 

KOM = (^OM^MORCI^ 1 ' 3 - 



We also examine Model B in the case where a larger 
complex, MOR,„CI c , represses Pr. Detailed formalism 
is given in the appendix. 

Figure [8] shows the plots for the extended Model B 
with (m, c) = (2, 1). This version of the model shows 
sharper transition between the MOR2CI regime and the 
CI2 regime for the form of CI protein as one can see 
in the lower graphs for CI ratio. As for the form of 
MOR, substantial fraction of MOR2CI appears only in 
the intermediate range of [CI] in contrast to the case 
of the MOR CI heterodimer in Fig|7] This is because 
there are not enough MOR's in the high [CI] region to 
form MOR2CI complex because it requires two MOR 
proteins. 

In the right two columns of Figj6] the parameters that 
give the Pr repression fold in the range [50, 200] are 
plotted versus the resulting Pl repression fold for Model 
B. The 2 x 10 5 parameter sets are chosen randomly over 
the range of [10 -7 , 1] for K l, ^om, K C i 2 , and #mor : ci- 
One can see that broad range of parameter sets satisfy 
the repression fold criterion for Pr, but the resulting re- 
pression folds for Pl are limited to the region larger than 
5,000 for (m,c) = (1, 1) and 50 for (2,1). 



3. Discussion 

3.1. Summary 

The bacteriophage TP901-1 has provided us with a 
conceptually new design of a genetic switch, in which 
the interaction between two antagonistic regulators, CI 
and MOR, is essential. Bistability between the immune 
and the ant-immune states has been demonstrated with 
a genetic switch that consists of the two divergently ori- 
ented promoters Pl and Pr, the two promoter-proximal 
genes, cl and mor, and only one of the three CI opera- 
tor sites Ol on a low copy number plasmid [8, 9]. The 
repression folds in the two states have been determined 
by in vivo measurements as around 1,000-fold for Pl re- 
pression in the immune state and around 100-fold for Pr 
repression in the anti-immune state. 

We constructed mathematical models for this cloned 
bistable system, assuming a putative operator Om to 
regulate Pr (Fig[TJ)). We assumed that Pl is repressed 
by CI2 bound to Ol whereas Pr is repressed by the 
MOR„,CI e :DNA complex on Om . We examined two 
types of models: one where MOR and CI interact only 
on DNA (Model A), and the other where MOR and CI 
form MOR m CI c complex in cytoplasm first and then the 



Parameter Ranges for Model A 



(m,c) (2,1) (2,2) 



K OL ~ 10- 2 10" 3 ~ 10- 2 

^om lO" 2 ~ 10-' ~ 10" 3 

Kg 2 > IP' 2 > IP' 2 



Table 1: Accepted ranges of parameters for Model A. The values are 
given in the unit of the concentrations that correspond to [CI] and 
[MOR] at the full activity of the promoter Pr and Pl, respectively. 



Parameter Ranges for Model B 



(m,c) (1,1) (2,1) 



K OL < 3 x 1(T 3 (3 x 1(T 6 ) ~ 3 x KT 2 

Kom < 5 x 10- 3 (10~ 4 ) ~ 3 x 10-2 

Kch — — 

^MOR:CI < 3 X 10 1 (10~ 3 ) ~ 1 



3.2. Validity of the Models 

In order to assess validity of the models, we have 
to determine our unit for CI concentration first. We 
employed the unit where [CI] is measured by the con- 
centration in the steady state with the full activity of 
Pl. This concentration should be compared with [CI] 
in the immune state of in vivo experiment on the sys- 
tem with around 10 copy-number plasmid containing 
the modified switch. We estimate this as follows by 
using the value of 300 nM for the CI concentration in 
the lysogenic/immune state of the wild lambda phage 
rf20ll ; First, we assume that this concentration is compa- 
rable with that for wild type TP901-1 with a single copy. 
Then, we multiply this by the following two factors: the 
factor 10 of the copy-number of plasmid, and the factor 
100 of the relative activity of Pr in our modified system 
in comparison with the wild type switch [9]. With these 
factors, we estimate that pRo, i.e., [CI] at the full Pr ac- 
tivity in the present system, could be well over 10 5 nM 
scale. 



Table 2: Accepted ranges of parameters for Model B. The values are 
given in the unit of the concentrations that correspond to [CI] and 
[MOR] at the full activity of the promoter Pr and Pl, respectively. The 
values in the parentheses are the lower limits when the fold criterion 
for Pl is restricted to 500 < pL(open)/ />L(closed) < 5000. The ranges 
for Kqi 2 cannot be set because the accepted values extends over the 
whole tested range. 



complex binds to Om (Model B). For each model, we 
tested bistability and performed parameter scans using 
the criterion that the repression folds should be consis- 
tent with the experiments. 

Our results are summarized as follows: For Model A, 

(i) the system shows bistability only when 2m — c > 2, 

(ii) the possible values for the operator affinities Kol 
and Kom are narrowly distributed, (iii) the possible 
dissociation constant Kch is much larger than the 
operator affinity A'ol due to the large repression fold 
for Pl, (iv) the possible value for Kch is bounded by the 
relatively large lower limit. The accepted ranges for the 
parameters are listed in Table Q] For Model B, (i) the 
bistability is robust due to the sequestration of CI by 
the MOR,„CI c complex formation, (ii) large parameter 
regions are allowed by the repression fold criterion, (iii) 
the Pl repression fold is bounded by the lower limit for 
the parameters that are consistent with the Pr repression 
fold: pL(open)/pL(closed) > 5000 for (m,c) = (2,1), 
and > 50 for (m, c) = (2, 2). The accepted ranges for 
the parameters are listed in Table [2] 



Model A. With this unit, we might be able to rule out 
Model A based upon the estimated values of A"ol and 
Kci 2 ■ The possible value of Kol in Model A is between 
around 10~ 3 and 5 x 10~ 2 (Table [TJ, but this contradicts 
the in vitro estimate of 28 nM for the CI concentration 
at which Ol is occupied by CI for 50% of the time[ 10]. 
The lower limit of A'ol ~ 10~ 3 in Model A should cor- 
respond to 100 nM or quite possibly even larger, but it 
is already well above 28 nM, i.e., the in vivo estimate 
for repressor-DNA affinity for TP901-1. 

We also found that the large repression fold of Pl en- 
tails K C h » ^ol for Model A. This means that CPs 
exist as monomers in cytoplasm and form Ch when 
they bind to Ol, but this is in contrast with many phage- 
encoded repressor proteins, such as those encoded by 
phage lambda, 434, and 186, which tend to exist as 
dimers or higher oligomers in solution Jill E2I EM l24ll : 
Actually, most of the 434 and lambda repressors exist in 
the dimeric conformation at nanomolar concentrations 



II23L 12111 . Our Model A challenges the presumption that 
the formation of dimers is a prerequisite for its specific 
DNA binding. 

Model B. We found broader distribution of parameter 
sets that satisfy the repression fold criterion for Model 
B. In particular, we did not find lower bounds for possi- 
ble Kch in contrast to the case of Model A. 

In the comparison of the two variants of Model B, 
our results show that the model with the formation 
of MOR2CI complex is more favorable than that with 
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MORCI. For the model with MORCI, the Pl repres- 
sion fold turned out to be always larger than 5,000 for 
the parameters that give the Pr repression fold between 
50 through 200. Such a high repression fold of Pl has 
never been actually observed in vivo. On the other hand, 
for the model with MOR2CI, the lower bound for the 
resulting Pl repression fold can be as low as 50, which 
covers the observed range of the Pl repression fold. As 
for the parameter ranges, this variant of Model B gives 
^om > 0.2 x 10~ 4 and /T M or:Ci > 10~ 3 when the re- 
pression fold for Pl is smaller than 5,000 (Fig|6]l. 



3.3. Experimental Test 

One of the distinguishing consequences of Model B 
is that the system with uncontrolled Pr can be bistable 
because of the sequestration of CI by MOR:CI complex 
formation. Even if Pr produces CI at a constant rate, 
there can be the two stable states: one with the repressed 
Pl and the other with the derepressed Pl. Such a mech- 
anism of bistability has been proposed by Francois and 
Hakim B12l 1 1311 as a theoretical possibility. Our study 
suggests this mechanism is employed in TP901-1. This 
may be tested experimentally for the genetic switch of 
phage TP901-1 by measuring the promoter activity of 
Pl in systems containing a functional mor gene and ex- 
pressing CI from uncontrolled Pr promoters at constant 
but various rates. Plotting the Pl activity of each sys- 
tem versus the uncontrolled Pr activity, one should find 
a characteristic feature for bistability as in Fig|9] where 
the Pl activity is doubled-valued for a certain range of 
the Pr activity. 



(m, c)=(2, 1), pL =l, P R =1 



3.4. Concluding Remark 

The genetic switching mechanism in TP901-1 is re- 
markably robust; The modified system studied here 
with only one operator Ol contains 100 times more CI 
molecules in its immune state than the wild type genetic 
switch with all of the three operators on plasmids, yet 
still shows bistability. Our model study suggests that the 
robustness of the genetic switch in TP901-1 is brought 
about by sequestration of CI through MOR: CI complex 
formation in cytoplasm. 
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K OL =7.69 xlCT 
^=4.18x10 



K a =1.55 xl0 J 




PR activity, or [CI] tota | 



Figure 9: The Pl activity (or [MOR] tC)la i) versus the Pr activity (or 
[CTJtotal) f° r the system with uncontrolled Pr. The solid lines show the 
stable states and the dashed line shows the unstable state. 



Appendix: Formalism for Model B with (m, c) = 

(2,1) 

In the appendix, we present some of the formulae we 
used for Model B with (m, c) = (2, 1). For this case, the 
promoter Pr activity is a function of MOR2CI concen- 
tration, 



P R([MOR],[Cl]) = P Ro 



and [MOR 2 CI] is given by 



1 



1 + [MOR 2 CI]/^om 



[MOR2CI] 



[MOR] 2 ■ [CI] 
iF 2 " ' 

"■MOR:CI 



(26) 



(27) 



with the dissociation constant A'MORiCi- 

The total concentrations of CI and MOR are now 

[OWa = [CI]+2[CI 2 ] + [MOR 2 CI], (28) 
[MOR] to tai = [MOR] +2[MOR 2 CI], (29) 

thus the corresponding equations with Eq.(l20b in the 
main text becomes 



[MOR] 



K 1 

"-MOR:CI 



4[crj 

^l + 8[CI][MORWXj 



MOR:CI 



(30) 
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Thus, both sides of the steady state condition Eq.© are 
now given by 



pR([MOR],[CI]) = pR 



1 +[MOR] 2 [CI]/£. 3 



-(31) 



OM 



[CI] to ,ai = [CI]+2^-+-([MOR] total -[MOR])(32) 

with [MOR] by Eq.® and [MOR] tota i by Eq.©. The 
effective affinity Kqm is defined as 



Kqm = (KomK mor . ci ) 



1/3 



(33) 
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The stability criterion for the steady state used in the manuscript is examined. Dynamical analysis 
shows that the criterion is valid as long as the [CI] to tai curve is a monotonically increasing function. In 
Model B, however, the [CI] tota i curve has a part with negative slope for some parameter region, in which 
case the stability cannot be determined only by comparing the slopes of the curves. The steady state at 
the left side could be unstable when it is loacated in the region where [CI] to tai is decreasing. 

1. Steady States 



The steady states satisfy the self-consistent conditions 

pL([CI]) = [MOR] tot ([MOR],[CI]) (1) 

pR([MOR],[CI]) = [Cr\ tot ([Cr\,[MORj), (2) 

with 

[CI] tot ([CI],[MOR]) = [CI] + 2[C/ 2 ] + [MOR • CI] (3) 

[MOR] tot ([MOR],[CI]) = [MOR] + [MOR CI]. (4) 



Here, [CI 2 ], [MOR ■ CI] are equilibrium concentrations of the protein complexes. 

We determined the steady solutions graphically by looking for intersections of the following 
two curves as a function of [CI], i.e. the production curve and the [CI] to tai curve that represents 
degradation rate: 

pR = pR^[MOR], [CI]j (5) 

[C7] tot = [Cr\ m \[Cr\,[MOR]j, (6) 

with [MOR] being a function of [MOR] tot and [CI] derived from the relation 
[MOR] tot ([MOR], [CI]), and [MOR] tot being given by pL([CI]\ 



[MOR] = [MOR]([MOR] tot , [CI]) = [MOR](pL([CI]), [CI]). 

1 



(7) 



2. Stability criterion 



In the text, the stability of the steady state is determined by the simple criterion. Let the slope of 
the production curve and the degradation curve be denoted by 

dpR d[C/] tot 
and 



d[CI] d[CI] 
respectively. Then the criterion is 

the steady state is stable if^<^gf. 

(8) 

dpR d[C/] tot 

the steady state is unstable if > 

y d[C7] d[C7] 

at the corresponding intersection. 

This criterion is simple and plausible, but based on the single variable picture although the 
system has at least two dynamical variables: [CI] and [MOR]. The full analysis for stability 
requires dynamical consideration. 

3. Dynamical Analysis of Stability 

The dynamics for the protein concentrations is given by the set of equations: 
j f [MORU = ^-(pL([CI])-[MOR] tot ) (9) 

^[C/] to t = —[pR(\MOR], [CI]) - [C7]tot)- (10) 

The total concentrations [C7] tot and [MOR] tot are given by eqs.© and ©. Basic assumption for 
this is that the equilibration among protein complexes in cytoplasm is much faster than the decay 
rates of CI and MOR: 1 /r c and 1 /r M . 

Consider the steady solution with [CI]* and [MOR]*, which satisfies eqs.® and ©. Suppose 
the steady state is perturbed by small fluctuation as 

[CI] = [CI]*+6[CI] (11) 
[MOR] = [MOR]* + 6[MOR], (12) 

and see if the small deviation will grow or decay in time. 

By inserting these into eqs.© and (flOl) . we obtain the equations for the time evolution of 6[CI] 
and 6[MOR], 
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T C \\dC) \dc))' Tc\\dM) \dM, 
where we employ the abbreviated notations: 

M t = [MOR] toU M=[MOR], C t = [C7] tot , C = [CI], L = pL, R 
We further abbreviate the notation as 

>dpL([CI])\* 



then eq. ([T3T) is expressed as 



f M fiC , M UM \ 



c t ,c j Q 



v 8M j 



d[CI] ) ' 
— (L c - M uc , Af fiAf 

— (^,c _ c t A, — (r,m — CV,m) 



etc. 



4. The criterion is always valid for Model A: 



In the case of Model A, [MOR ■ CI] = 0, thus we have 
Q M = 0, M uM = l, M uC = 0, 
then eq. ([T4l) becomes 



' 0, 1 ^ 




(*?)=' 




V 



i(/? >C - C f , C ), M , 



8M 



Now we assume the solution as 
SC, 8M oc e wt 
then we have 



6C 
6M 



0. 



The condition that this equation has non-zero solution gives 



- Q c of + 



—(Rc - C t A - —Ctjc to + \l c r m + R,c - C t A = 

3 



If all the solutions co have a negative real part, the state is stable, whereas the state is unstable if 
there is a solution with a positive real part. 



Note that the slopes of the production curve and the degradation curve are given by 



dpR _ ( dpR([MOR],[CI]) \* 

= I I — L c A M + ( 

dCI \ d[C7] / ' ' 



dC/tot ^ ( d[CI] tot ([CI]) \* 
dCI \ d[CI] j M 

Using these, eq. ([T71) becomes 



(18) 
(19) 
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Note that 

L,c, R,c> R,m < 0. 

(1 ) The case ofr M <^ t c . One solution is of order 1 /t m and the other is of order 1 /t c . 
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This shows the criterion © is valid. 

(2) The case ofr c «: t m - One solution is of order l/r M and the other is of order l/r c . 
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The criterion ([8]) is also valid because 
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(3) General case. The criterion © can be shown to be valid because eq.dTTT) always has two real 
solutions, co i and co 2 , and the sum of the two solutions is negative when the slope of the production 
curve is less steep than that of the degradation curve, i.e. 

dpR dC/tot „ 
coi + co 2 < when — — - — — - < 0. 

dCI dCI 



The positivity of the discriminant D of eq.U7\ 

n 2 



D = 



1 1 

— {R,c - Q,c\ Qc 

+4Q C — [L C R M + R c - C uc ] 
t mTc 

% c _(i_±)c,; 2 



thus eq.([T71) has two real solutions. 
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5. The criterion is valid for Model B as long as the degradation curve has pisitive slope, but the state may be 
unstable otherwise. 



In this model, the existence of MOR ■ CI makes the expressions for the slopes of the production 
curve and the degradation curve a bit more complicated: 
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where the variables that kept constant upon partial differentiation are explicitly indicated as 
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whenever it could be ambiguous. In the derivation, we have used the relations 
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5.1. Stability of the system with PL with externally controlled [C/] to t 

First, we will examine the stability of the system without PR, CI being provided externally. The 
system is shown to be bistable for some parameter region. 

Based upon the approximation that the relaxation in the solution is much faster than the protein 
production rate by PL, we consider the system where [CI] and [MOR] satisfy 

[CI] tot = [CI] tot ([CI],[MOR]) (23) 

j f [MORU = ^-(pL([CI])-[MORU), (24) 

thus the deviation from the steady state follows 
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with the notation 

A = C uC M tM - C tM M uC > 0, 

whose inequality can be shown from the actual expressions of [CI] tot and [MOR] tot , © and ©. 
Therefore, we have 
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5.2. Stability of genetic switch with PL and PR: 

The growth rate co is determined by the characteristic equation 
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which can be expanded as 
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This can be put in the form 
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This is almost the same with the corresponding equation for Model A eq.(l20l). 
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(1) The case o/tm <^ Tc- One solution is of order 1 /t m and the other of order 1 /r c : 
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Therefore, the criterion © is valid as long as the slope of degradation curve is positive, i.e. 
(dC/tot/dC/) > 0, but the state is always unstable for the negative slope for the degradation curve, 
i.e. (dC/tot/dCT) < 0. 



(2) The case ofr c «; r M . 
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In this case, the criterion © is always valid because 
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(3) General case. We can show eq.(l30l) always has two real solutions, oj\ and a> 2 , and the steady 
state stability is determined from the sign of u)\ + co 2 : 
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This means that the stability cannot be determined only from the slopes of the production and 
the degradation curves in the case the degradation curve has a negative slope. 

Appendix: The positivity of the discriminant D of eq. 
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Thus, D is 
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The second line of the last expression is shown to be positive as follows: 
Note that 

Ctjc = 1 + C 2 ,c + MC, C , C tM = MC M , M tjC = MC, C , M tM = 1 + MC, C , 
R, c = R' ■ MCjc, R, M = R' ■ MC M , 

where the abbreviated notations are used, 
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Therefore, 
D > 0. 
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